Entangled light in transition through the generation threshold 
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We investigate continuous variable entangling resources on the base of two-mode squeezing for 
all operational regimes of a nondegenerate optical parametric oscillator with allowance for quantum 
noise of arbitrary level. The results for the quadrature variances of a pair of generated modes are 
obtained by using the exact steady-state solution of Fokker-Planck equation for the complex P- 
quasiprobability distribution function. We find a simple expression for the squeezed variances in the 
near-threshold range and conclude that the maximal two-mode squeezing reaches 50% relative to 
the level of vacuum fluctuations and is achieved at the pump field intensity close to the generation 
threshold. The distinction between the degree of two-mode squeezing for monostable and bistable 
operational regimes is cleared up. 

PACS numbers: 03.67.Mn, 42.50.Dv, 42.50.-p 



Recent years have witnessed an increasing interest in 
quantum information theory dealing with entanglement 
of continuous variable (CV) systems CV entangled 
states of light hold the key for quantum communications 
at the high-intensity level, as CV analogues of various 
protocols developed originally in the framework of dis- 
crete quantum variables have been established (see, for 
example Q). An important example provides en- 
tangled Einstein-Podolski-Rosen (EPR) state of orthog- 
onally polarized but frequency degenerate beams, effi- 
ciently generated via nonlinear optical process of para- 
metric down-conversion. It has been pointed out in Q 
and has been demonstrated experimentally in Q for CV, 
employing sub-threshold nondegenarate optical paramet- 
ric oscillator (NOPO). Then a CV entanglement source 
was built from two independent quadrature squeezed 
beams combined on a beam-splitter The generation 
of CV EPR entanglement usingan optical fibre inter- 
ferometer is also demonstrated It should be noted 
that so far most experimental realizations of the CV en- 
tanglement on NOPO's have only been operated below 
threshold. A natural next step is the extension of these 
investigations to laser-like systems generating entangled 
bright-light states. Unfortunately, both theoretical and 
experimental studies of these problems are very com- 
plicated and only rare examples are known up to now. 
Many efforts have been devoted to the study of inten- 
sity correlated twin beams from NOPO above threshold 
0. The conditional generation of sub-Poissonian light 
from bright twin beams in the CV regime is experimen- 
tally demonstrated in fsj. Further experimental studies 
of bright two-mode entangled state from cw nondegener- 
ate optical parametric amplifier have been made in Refs. 
0. The generation of CV polarization entanglement by 
mixing a pair of polarization squeezed beams is realized 
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in Tol. 

The purpose of this paper is to investigate physi- 
cal properties and the presence of CV entanglement for 
NOPO in its transition through the generation threshold 
as well as in the regime of lasing. One of the princi- 
pal problems in this study is the description of quantum 
fluctuations. In most theoretical works nonclassical ef- 
fects and entanglement resources of nonlinear quantum 
systems are usually described within linear treatment of 
quantum noise. It is obvious that such approach does not 
describe the critical ranges (threshold, point of multista- 
bility, etc) where the level of quantum noise increases 
substantially. We use a more adequate approach within 
the framework of exact nonlinear treatment of quantum 
fluctuations via the solution of the Fokker-Planck equa- 
tion for the quasiprobability distribution function. De- 
riving the quasiprobability functions so far has been per- 
formed only for a few simple models (see and ref- 
erences therein). For the NOPO the exact steady-state 
solution of the Fokker-Planck equation in the complex 
P-representation first was obtained in [121. We will use 
the generalized form of such solution [iSj which involves 
also dctunings of the modes and hence admits bistability. 

There are various questions that emerge in the study 
of these problems. Because NOPO displays both monos- 
table and bistable regimes it is important to understand 
how the properties of entanglement depend on the opera- 
tional regimes. What are peculiarities of entanglement in 
the vicinity of threshold? Will entanglement take place 
in the regime of lasing or how far it can be extended 
into the high intensity domain? Answering these ques- 
tions is extremely important for a deeper understanding 
of quantum entanglement, and also from the perspective 
of creating an entangled light laser. 

We consider a type-II phase matched NOPO with 
triply resonant cavity that supports the pump mode at 
the frequency co^ and two orthogonally polarized modes 
at the same frequency lui ^ lu2 ^ ^3/2. The pump mode 
is driven by coherent field at the frequency lul — i-^s- The 
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relevant interaction Hamiltonian is 

H = ihn (oj'^a^as — 01020;^) + (1) 

3 

2 = 1 

where (i = 1, 2, 3) are the operators of the modes uji, k 
is the couphng constant, E = \E\ e'*^ is the complex am- 
plitude of the pump field. The last term in |^ describes 
mode damping in the cavity in terms of the reservoir op- 
erators Ti and Tf , which determine the damping factors 
7i. We take into account the detunings of the cavity 
A3 = lol — i^3, and Ai = A2 = A = i-ul/2 — uji^2- We 
also assume perfect symmetry between the orthogonally 
polarized modes provided that they decay at the same 
rate: 71 = 72 = 7- 

It is well known that linear treatment of NOPO is not 
self-consistent due to the phase diffusion. According to 
this effect the difference between the phases of the sig- 
nal and idler modes, as well as each of the phases can 
not be defined in the above threshold regime of gen- 
eration. On the whole, the well founded linearization 
procedure cannot be applied for this system. Neverthe- 
less, the linearization procedure and analysis of quantum 



fluctuations for NOPO become possible due to the addi- 
tional assumptions about temporal behavior of the dif- 
ference between the phases of the generated modes [T^ . 
Within this assumption the two-mode squeezing spec- 
tra are obtained for both below- and above-threshold 
regimes. These results are inapplicable in the near- 
threshold as well as in the bistable ranges of generation, 
where a more accurate treatment of quantum noise is 
necessary. Our analysis of NOPO is based on the steady- 
state exact solution of the Fokker-Planck equation. It 
is obvious, that in the framework of this approach we 
avoid both difficulties of the linear theory. The price one 
has to pay for this advantage is the impossibility to per- 
form temporal description of quantum fluctuations and 
hence squeezing spectra cannot be calculated. Neverthe- 
less, our analysis does not require spectral description of 
squeezing and therefore we do not encounter any princi- 
pal difficulties. 

We start with the Fokker-Planck equation of the sys- 
tem in the so called complex P-representation of the den- 
sity matrix and in the case of high cavity losses for the 
pump mode (73 ^ 7), however, in the operational regime 
the pump depletion effects are involved: 
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P(a, t). 
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Here a = (ai, aj^, a2, aj), at, af [i = 1,2) are the in- 
dependent complex variables corresponding to the oper- 
ators Oi, a+, and 7^- = 7^ - lAj, (7^ = 72 = 7)- 

The normally-ordered moments of time-dependent op- 
erators are calculated through the P-quasiprobability dis- 
tribution function as 

{at{tfa,{t)'at{tra2{tr) = (3) 

= j da^daida2da2P{a, t)a'^'' a[ a 2 a2 ■ 



Below we consider only the stationary steady-state 
regime and drop the time dependence of operators when 
calculating one-time expectation values. Using the 
steady-state solution of the equation ^ , obtained in [l3l . 
and the method of integration in the complex plane 
we find: 
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{a+"'a'['a+"a^) = (a+"a>+™a^) , if n > to, 



(6) 



where s = E/k, p = |2e|^ is the scaled pump intensity, 
A = {x)j := x{x+l)...{x + j -1), {x)o = 1, and 



j=o (A + l)_ 



3=0 



(7) 



In addition to these expressions we note that the system 
under consideration has the foUowing property, condi- 
tioned by its most general symmetries: 

I) =0, iik-l^m-n. (8) 
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This property is the consequence of the rotational sym- 
metry of the Hamiltonian ||IJ U (9) HIJ-^ {9) = H for 
any 9, where U (9) = exp [i9 [afai — 0^02)] . Since the 
Lindblad part of the master equation is invariant with 
respect to such transformation too, then the steady state 
density operator of NOPO must commute with U (9): 
U (9) pU~^ (9) — p, whereupon we arrive at the relation 
®. 

We will use only the photon number and the (0102) 
moment in the further discussion, so we write their ex- 
pressions explicitly: 



2^ ■ j 

= (a+ai) = (0+02) = ^ E 1 — 2 

j=i (A+1) 



(0102) = 
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2^VP,tt |(A + 1)^. 
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We note that the state generated in NOPO is non- 
Gaussian one, i.e. its Wigner function is non-Gaussian 
|15|. So far, the inseparability problem for bipartite non- 
Gaussian states is far from being understood. On the 
theoretical side, the necessary and sufficient conditions 
for the separability of bipartite CV systems have been 
fully developed only for Gaussian states |ljj . To charac- 
terize CV entanglement we have chosen the inseparabil- 
ity criterion based on the total variance of a pair of EPR 
type operators. For a pair of optical beams generated in 
NOPO this criterion characterizes the entanglement in 
terms of quadrature operators 



Yk - Yk{9k,t)^Xki9k~^,t), 
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(fc = 1,2) and due to the mentioned symmetries is re- 
duced to the following form: 

V{9i,92) := V (Xi -X2) = V (Fi + Y^) < 1. (12) 

Here V{x) is a convenient short-hand notation for the 
variance (a;^) — (x) . Inequalities (|12|) require the vari- 
ances of both conjugate variables Xi — X2 and Yi + Y2 



to drop below the level of vacuum fluctuations. Since the 
states of NOPO are non-Gaussian, the criterion H12() is 
only sufficient for inseparability. 

At this point we must note the difference between the 
focus of our paper and most of the preceding work de- 
voted to the study of two-mode squeezing. It is an estab- 
lished standard to describe squeezing with the spectra of 
quantum fluctuations of the considered variables. The 
recent experiment on spectral investigation of criteria for 
CV entanglement was presented in the paper ■ Unlike 
that, we aim to analyze the separability properties of the 
system solely with the help of the criterion H12|l , which in- 
volves a single integral characteristic of the NOPO quan- 
tum state. For the completeness of the results it is use- 
ful to remind how intracavity variances are related to 
the spectra of squeezing of output fields. For this, let 
us consider the special scheme of generation, when the 
couplings of in- and out-fields occur at only one of the 
ring-cavity mirrors. For the case when only the funda- 
mental mode is coherently driven by the pump field, while 
subharmonic modes are initially in the vacuum state, we 
have for the output fields a°"\t) = ,y2^a,{t), (i = 1,2). 
The spectra of two-mode quadrature amplitude squeez- 
ing is 



SiLO,9i,92) 
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t + T))dT. 



Here (0i , 02 , = [^1(^1, " ^2(^2, i)], and 
t is sufficiently large {t ^ 7"^) to ensure steady state, so 
that S{uj,9i,92) does not depend on t . 

The variance V {9 1,92) is expressed with the above 
spectra in the following way: 

271/(01,02) = 7^ / SiLU,9i,92)dLu. (14) 

Let us turn to the integral variances. With an ap- 
propriate selection of the phases 9k of the quadrature 
amplitudes, namely 

01 + 02 = arg (0102) , (15) 

we arrive at the following minimum value of ^(01, 02) 

V^nin = l + 2(n- |(aia2)|). (16) 

Detailed physical analysis of the variance IjlGI) on the 
base of the exact expressions ©, UlUf) seems to be a 
rather complicated task. To proceed further we trans- 
form the formula (|10|l to a more appropriate form 



^a,a2)^e--^(A*+2n+'Pp]. (17) 
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FIG. 1: Scaled mean photon number versus dimensionless 
amplitude of the pump field Es = (2K/773) |_E| for monos- 
table (curve 1, A/7 = 1), interjacent (curve 2, A/7 — 3) and 
bistable (curve 3, A/7 = 7) regimes {n/y — 0.5, 73/7 = 18). 
Dot curve visualizes the classical analysis for the bistable case. 
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FIG. 2: Minimized variance V {Xi — X2) versus dimension- 
less amplitude of the pump field Es for the same parameters 
as for Fig.Q 

What is important is that (0102) is expressed with n 
in a simple enough form, much easier to understand than 
its expansion in powers of p. The advantage of that form 
is due to the intuitiveness of the behavior of n in vari- 
ous feasible operational regimes. Moreover, semiclassical 
solution for the mean photon number, where applicable, 
can be employed to obtain a semiclassical solution for 
(0102). To this end, we remind relevant results of the 
semiclassical approximation. In the regime below thresh- 
old p < p^^ the excitations of modes are at the level of 
spontaneous noise, and above threshold they have the 
form 

Uci - V- Re A + ^/p- (ImA)'j . (18) 




E 
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FIG. 3: Illustration of the critical region squeezing for small 
values of k/7 ratio. K/7 = 0.5 for the solid line, and K/7 = 
10"'' for the dot fine. 



The case ReA > corresponds to monostable dy- 
namics with threshold = |A|^ (The corre- 
sponding threshold value of the pump field = 
1(7 — iA) (73 — iAs)! /k.) The opposite case ReA < 
is bistable with threshold p*j^ — (ImA)^ and the stabil- 
ity region of the zeroth solution extending up to p = 
|A|^. For completeness, we present in Fig. ^ quantum- 
mechanical results for mean photon number of the modes 
@ for three values of detunings corresponding to monos- 
table, bistable and interjacent (ReA = 0) regimes. For 
the bistable dynamics the critical region, i.e. the range 
of the pump intensity, where the system passes onto the 
nontrivial classical branch, lies between p^'* and |A|^. The 
quantum critical region of monostable dynamics is in the 
vicinity of the generation threshold = |A|^. 

We are now in a position to study the entanglement 
effects and will state the main results of the paper. 
What is important is that we find the quantities l|^- l(Tn|) 
through the exact steady-state solution of the Fokker- 
Planck equation and carry out an exact quantum statis- 
tical analysis. 

Fig. 121 illustrates the dependence of Vrain on the pump 
amplitude. One can see that for monostable dynamics 
(curve 1) entanglement is realized in the entire range of 
pump intensities. In all cases maximal degree of two- 
mode squeezing Knin = 0.5 is achieved within the critical 
region. Above the critical point mean photon numbers of 
the modes increase considerably and variance T4iin starts 
to increase too. For the bistable case (curve 3) growth 
of V^iin is much faster and larger so that the sufficient 
criterion for inseparability \i2\ is not fulfilled. Note that 
T^nin has a sharp peak in the critical range of bistabil- 
ity. Anyway in the far-above threshold region the CV 
entanglement of the system is guaranteed. More exactly, 
the asymptotic value of Knin for all detunings is equal to 
0.75. 

Attentive reader of this article must be somewhat con- 
fused by the unnatural value of K/7 ratio selected to plot 
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the figures. We emphasize that the choice of parame- 
ters is conditioned merely by illustrative purposes. For 
smaller, experimentally available values of the parameter 
k/7 ~ 10~^ 10~^, the behavior of Kiiin is not changed 
qualitatively in comparison with the results of Fig. |21 
To illustrate this fact we plot in Fig. 13 the behavior of 
Vmin in the interjacent regime and for two values of K/7. 
One can see that for the more realistic parameter value 
k/j = 10~® behavior in the vicinity of the critical point is 
more abrupt, which is just a result of critical region being 
much narrower. For the bistable dynamics this narrow- 
ing of the critical region makes the peak of Knin (Fig. 
12 curve 3) much sharper, the rest of the plot remaining 
unchanged. 

It is remarkable that equations (|16|l . (|17|l allow to ex- 
plain qualitatively and even quantitatively the numerical 
results obtained. 

1. In the far below-threshold regime mean pho- 
ton number can be represented linear in p: n = 
p/2|A+l|^ <C 1. Using this formula it is easy to check 
that the condition Vmin < 1 is fulfilled for any A within 
the entire domain where the weak pump approximation 
is valid. 

2. In the far above-threshold range {p ^ IM'^) '^^ '^^^ 
make use of the semiclassical expression (I18f) together 
with H10|l to study the variance T^nin- To this end we 
write n — nd+Sn, and neglect Sn only where that doesn't 
lead to loss of accuracy in the formula After some 
algebra (see Appendix^) we arrive at the following ex- 
pression: 

lim Klin = 0.5 - 2Sn. (19) 

p — ^00 

Straightforward, but complicated analytical calculations 
on the formula © show that 5n —0.125 in the 
limit p ~f 00, which leads to the asymptotic value 
Vmin = 0.75 < 1 (find in Appendix 151 several hints re- 
garding our analysis of 6n). Therefore, as our analysis 
shows with allowance for quantum fluctuations of arbi- 
trary level, CV entanglement is always achieved in NOPO 
above threshold. What is remarkable, is that a very small 
quantum correction 6n to the semiclassical intracavity 
photon number plays an essential role to the production 
of entanglement in high-intensity level. Obtained result 
seems interesting as provides the example of preserving 
CV entanglement in the high-intensity domain. 

3. In the near-threshold and bistability regimes an- 
alytical study is too complicated. We have carried out 
strict analysis only for the monostable case (AA3 < 773), 
assuming k/j <C 1. We perform expansion of Vmin 
around the minimum point Pmin in powers of parame- 
ter s := |A| ' (s ^ keeping only the significant 
terms (technics similar to those described in Appendix IbI 
is used). The result, expressed through the pump field 
intensity / ~ , is the following: 

K,i„ = 0.5 + cVi(c). + ^ ' , (20) 

s \ Ith J 



where 

/min = Ith [1 + h (C) S] , (21) 

c = v/cos(argA), /i(c) > 0.0164, /2(c) = 0.113 -I- 
0.00221c -0.330c2 + 0.371c3-0.132c^ /3(c) = -2.219 + 
0.217c + 2.83c'^. Note that for fixed damping rates 7^ 
c is a decreasing function of detuning, and c = 1 at 
exact resonance. Formula H20|) is valid for c^ ^ s 
(i.e. not too close to the interjacent regime c = 0) and 

\{I-Imin)/Ith\ « 

Let us discuss the results (jJOl, in more detail. 

Positiveness of /i(c) means that the minimal value of 
Klin doesn't drop below 0.5. This result differs from the 
studied case of the perfect two-mode squeezing generated 
in an undamped NOPA, where Knin vanishes. Evidently, 
in our system the larger limiting value for V,nin is due to 
dissipation and cavity feedback effects. Then we see that 
the point of maximal two-mode squeezing Jmin is located 
close to the generation threshold Ith- Function /3(c) is 
negative for values of c close to zero, however it increases 
and becomes positive as c approaches 1. Thus the point 
of maximal two-mode squeezing can be located both be- 
low and above the generation threshold of monostable 
dynamics. 

Formulae (|20|l and (|21fl are an example of a result that 
couldn't be obtained with the help of linear treatment 
of quantum fluctuations. Although the linearized the- 
ory gives correct expressions for Vmm in the below- and 
above-threshold ranges, it not only fails to link those solu- 
tions through the narrow range of the generation thresh- 
old of monostable regime, but also turns entirely invalid 
in the interjacent and bistable regimes. 

In the light of experimental investigations of CV 
quadrature entanglement in NOPA and NOPO [l8l|. 
the obtained value for the highest achievable squeezing 
of 3dB (Klin = 0.5) may seem inadequate. Indeed, the 
measurements of the frequency spectra of the variances 
of the output quadrature-phase amplitudes reveal up to 
3.7 dB squeezing for NOPA 9], and up to 4.9 dB squeez- 
ing for NOPO operating above threshold 18]. However, 
we remind that our results pertain to the full squeezing 
of the system and not the spectral component squeez- 
ing. Moreover, as noted in |l3, states with even perfect 
squeezing of a certain spectral component may be non- 
squeezed in the full sense. Note, that complete study of 
that question is beyond the scope of this paper. 

Our results allow us to explain qualitatively the growth 
of Klin in the critical region of bistability. We remind 
that the minimization of the variance V is performed 
through a specific selection of phases of the quadra- 
ture amplitudes Xk and Yk - see formula (|15|) . In the 
bistable regime the state can be regarded as a mixture 
P = PwPw+PsPs of two states corresponding to the weakly 
and strongly excited semiclassical stable solutions. What 
if both pw and ps are highly squeezed two-mode states, 
but the minimizing phases are rather different for them? 
In this case the minimal variance for their mixture will be 
achieved for an intermediate value of the phases, which 
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FIG. 4: O = arg(aia2) — <1?b versus dimensionless amplitude 
of the pump field Ea for the same parameters as for Fig. Q 



is good for none of them. The analysis of the phase of 
(0102) suggests that this is indeed the reason. Fig. ^ 
shows that the minimizing phase changes very abruptly 
in the bistable critical region. In this context we note that 
in the critical bistable regime the state is extremely differ- 
ent from Gaussian and therefore the condition V^nin > 1 
does not witness against entanglement since the criterion 
Elprovides only sufficient condition for CV entanglement 
in this case. 

In summary, we have investigated CV entangling re- 
sources for all operational regimes of nondegenerate OPO 
in the most basic and explicit way through the P-complex 
probability distribution. We have studied the entangle- 
ment as two-mode squeezing and have shown that entan- 
glement is present in the system for quantum noise of 
arbitrary level and, what is remarkable, for a wide range 
of intensity of the coherent driving field including the 
far-above-threshold limit. The numerical and analytical 
analysis for the quantum critical range of NOPO has pro- 
vided a clear distinction between the degree of two-mode 
squeezing for monostable and bistable regimes. Most fa- 
vorable is the monostable regime where CV entanglement 
can be maintained for all values of the pump intensity and 
for experimentally available parameters. In our study we 
have not analyzed any entanglement measure which is ex- 
tremely difficult to handle analytically for non-Gaussian 
states |23|- This topic is currently being explored and 
will be the subject of forthcoming work. We hope the re- 
sults obtained can serve as a guide for further theoretical 
and experimental studies of bright CV entangled light. 
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APPENDIX A: ASYMPTOTIC FORM OF Knin 

Supposing that for p 3> |A|^ Sn is monotonous and 
5n <^ Uci (these assumptions are justified during a later 
analysis) it can be easily shown that n w ^y/p and 



p dn 1 
lim — — — — , 

p^oo n op 2 



(Al) 



The following exact expression is also not difficult to 
obtain using formula ((T7|l : 



\{a,a,)\=nJl + ^(2Sn+Pp 
y p \ n op 



where 



Q = 2Jp- {ImAy + 25n + 



p dn 
n dp 



(A2) 



(A3) 



Q ~ '2'y/p, when p ^ 00. Therefore ljA2|l can be ex- 
panded into series, whereupon wc find 

\{a,a,)\ = n + ^ (2^^ + + 0(p-^), (A4) 

This immediately yields the asymptotic behavior H19|) 

of Klin ■ 



APPENDIX B: ASYMPTOTIC VALUE OF 5n 

We decompose the value of 2nci into the integer and 
fraction parts: 2nci = /i -I- ^, with fi — E [2nci\ and ^ = 
F{2nci}- After that we can write: 



nci 



F 
2N 



2' 



where 



(Bl) 



(B2) 



3=0 



and N and Nj are defined with formula Q. 

It can be shown that A^^ is the largest term in N and 
that the values of N and F can be calculated with suf- 
ficient precision by summing only about 0(^//I) terms 
around Nf^. Within that range of indices the term A'^i+j 
can be represented as below: 



~ exp 



(B3) 



Summation of such terms can be safely replaced with 
integration, with bounds extended to infinity. We finally 
obtain the following expression for N: 



N 



/TT/i exp 



4a* 



(B4) 
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Calculation of F requires more accurate treatment, 
since it contains positive and negative parts. We en- 
sure against errors due to subtraction of two large values 
rewriting F as follows: 

o(v7J) 

P= jiN^+,-N^^,). (B5) 



The difference N^+j — N^-j can be simplified to an 
appropriate approximate form, allowing to replace the 
sum with integration. Proceeding like with N we arrive 
at6n^ - 0.125. 

p — >QO 
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